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Abstract 

The rates or formation and concentration distributions of a dimer re- 
action showing hysteresis behavior are examined in an ab initio chemical 
reaction designed as elementary and where the hysteresis structure pre- 
cludes the formation of transition states (TS) with pre-equilibrium and 
internal sub-reactions. It was discovered that the the reactivity coeffi- 
cients, defined as a measure of departure from the zero density rate con- 
stant for the forward and backward steps had a ratio that was equal to the 
activity coefficient ratio for the product and reactant species. This sur- 
prising result, never formally noticed nor incorporated in elementary rate 
expressions over approximately one and a half centuries of quantitative 
chemical kinetics measurement and calculation is accepted axiomatically 
and leads to an outline of a theory for the form of the rate constant, in 
any one given substrate - here the vacuum state. A major deduction is 
that the long-standing definition of the rate constant used for over a cen- 
tury and a half for elementary reactions is not complete, where previous 
works almost always implicitly refer to the zero density limit for strictly 
irreducible elementary reactions without any any attending concatenation 
of side-reactions. This is shown directly from MD simulation,where for 
specially designed elementary reactions without any transition states, den- 
sity dependence of reactants and products always feature, in contrast to 
current practice. It is argued that the rate constant expression without re- 
actant and product dependence is due to historical conventions for strictly 
elementary reactions. From the above observations, a theory is developed 
with the aid of some proven elementary theorems in thermodynamics, and 
expressions under different state conditions are derived whereby a feasible 
experimental and computational method for determining the activity coef- 
ficients from the rate constants may be obtained under various approxima- 
tions and conditions. Elementary relations for subspecies equilibria and 
its relation to the bulk activity coefficient are discussed. From one choice 
of reaction conditions, estimates of activity coefficients are given which 
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are in at least semi-quantitative agreement with the data for non-reacting 
Lennard- Jones (LJ) particles for the atomic component. The theory de- 
veloped is applied to ionic reactions where the standard Bronsted-Bjerrum 
rate equation and exceptions to this are rationalized, and by viewing ion 
association as a n-meric process, ion activity coefficients may in principle 
be determined under varying applicable assumptions. 

Keywords: [1] elementary reaction rate constant , [2] activity and reactivity 
coefficients, [3] elementary and ionic reactions without pre-equilibrium. 
AMS Classification: 80A10, 80A30, 81T80, 82B05, 92C45, 92E10, 92E20. 

1 Introduction 

Previous work has detailed a hysteresis model of a simple dimer reaction ^JEDE! 
where the coordinates for molecular formation rj and breakdown rb are not at 
the same vicinity, as shown in Fig.QJ. If both these points coincided, then one 
could in principle define a volume region about the point rj = that would 
serve as the transition state pre-equilibrium TS, suggesting a composite reaction 
such as 

2A^TS^A 2 (1) 

which is therefore not strictly elementary since JQ) is a summation of elementary 
steps to yield a net reaction. The current model precludes such a possibility 
, implying a strictly elementary reaction process. Previous models, especially 
that of Eyring postulated the pre-equilibrium TS which was also interpreted 
experimentally 0] but the most recent developments in some cases by-passes 
this development Ej • In the theory that follows, the TS refers to the point 
in state space whose neighborhood always includes product and reactant states. 
The method for treating changes of potential via switch mechanisms were also 
outlined together with an algorithm which was implemented to conserve 

energy and momentum at high energy regions with steep potential gradients. 
This work focuses primarily on the form of the elementary rate constants for 
the reaction, and presents a way of determining activity and reactivity coeffi- 
cients; the two concepts are different, and will be denned in this work, where 
under certain conditions, the activity coefficients might be determined. These 
two types of coefficients are discriminated by first introducing an ansatz of the 
dependency of the reactivity coefficients on the rate constants and comparing 
the results with the actual kinetics of the simulation, where it was discovered 
that there were variations of the rate constant with concentration. Section © 
and its subsections give details of the methods used to determine the extent 
of variation. The ratio of the reactivity coefficient ^ agreed numerically with 
the results from the activity factor $ e , which is determined independently by 
extrapolation of the variation of the concentration equilibrium ratio determined 
directly from simulation. This observation is used to develop a theory for the 
form elementary reaction rate constants take, which does not accord with stan- 
dard assumptions. The above objectives are realized by first describing the 
model of the elementary reaction in Section The empirical result of the 
equivalence of the reactivity and activity coefficient ratios is developed in Sec- 
tion |PJ and the subsections where the consequences are discussed in depth. 
The uniqueness of the activity coefficient and various other associated results 
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appear in Section (|3.3l) . Verification of the theories concerning the role of ac- 
tivity coefficients in elementary reactions is given in reference to results from 
the literature for non-reacting LJ particles, where it is argued that the residual 
Helmholtz energy A res is a better measure for computing the activity coefficient 
for multi-component system than the residual free energy G res (Sections ()3.4.U 
^ . In either case, the activity coefficient for the unreacted particle increases 
with system concentration and is positive. Of two possibilities for the estimate 
of the activity coefficient derived from energy considerations of the dimer and 
single particle trajectory along the reaction coordinate, only one accords with 
observation derived from simulation data from the literature for the activity 
coefficient of the atom, with semi-quantitative agreement. From the form of the 
activity/availability coefficients, rudimentary mechanisms involving single and 
double stages (Sections I|3.6.1l3.6.3l) ')arc proposed for elementary reactions. It 
is suggested that the 2-stage mechanisms is the more appropriate, and an inter- 
pretation of the Bronsted-Bjerrum (BB) rate equation is made on the basis of 
the two-stage model (Sections (13.6. 413. 6.5j) y The conclusion (Section(@})brings 
into focus all the details of the preceding sections. 



2 The Model 

The simulation model is a dimeric particle reaction 

2A ^ A 2 (2) 

at the Lennard- Jones (LJ) supercritical regime (T* = 8.0, 0.03 < p < 1.1) 
in a range of equilibrium fluid states. Details of the mechanism have already 
been described and will not be repeated here 0|H]- In the current study, the 
potentials given in Fig. ^ are used. In this model, 2 free atoms '"react"' at rt 
where a switch mechanism converts the potential to the harmonic intermolecular 
potential. The molecule '"breaks'" at where the potential reverts to the LJ 
type. At these points, a specialized algorithm was used to conserve energy and 
momentum and this too has been described in detail Ej ■ The free atoms A 
interact with all other particles (whether A or A2) via a Lennard- Jones spline 
potential and this type of potential has been described in much detail [§| . An 
atom at a distance r to another particle possesses a mutual potential energy 
ulj where 

ulj = 4e (-) ^ - (-) 6 for r < r s (3) 



u LJ = a i0 {r - r c ) 2 + b tj {r - r c f for r s < r < r c 

ulj = for r > r c 

and where r s = (26/7) At 77, the potential is switched to the molecular 

potential given by 

u(r) = u vib (r)s(r) + u LJ [1 - s(r)] (4) 

where u V ib(r) is the vibrational potential given by eq.JHI below and the switching 
function s(r) has the form given by eq.JSJ). 

s{r) = j- — ttt (5) 

1 + 
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Figure 1: Potentials used for this work 



s(r) 
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if r < r sw 
for r > r„ 



. The switching function becomes effective when the distance between the atoms 
approach the value r sw (see Fig. (|Tp . The intramolecular vibrational potential 
Mt,ih(r) for a molecule is given by 



u vi b{r) = u + -k(r - r ) 2 . 



(6) 



LJ reduced units are used throughout this work unless stated otherwise. The 
details of the parameters have been given elsewhere |7] but they include the 
following: 

uq — —10, ro = 1.0, k ~ 2446 (exact value is determined by the other input 
parameters), n — 100, r/ = 0.85, = 1.20, and r sw = 1.11. 



3 Thermodynamic results from equilibrium mix- 
tures 

Details of the runs have been described elsewhere 0. Typical runs of 10 mil- 
lion (10M) time steps were performed at each general system particle density 
p (where p refers to the particle which is either free or part of a molecule), 
where the first 200,000 steps were discarded so that proper equilibration could 
be achieved for our data samples. The sampling methods have been previously 
described [5] where sampling of all data variables were done each 20 th time step 
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and where the data were averaged and written into a dump file of 100 dumps 
for the 10M million time steps. The averaged values in each dump were further 
averaged to yield the final averages and standard errors. Dynamical quantities 
however had to be sampled at each time step St* = 0.00005. In view of the ab- 
normally high temperatures -not hitherto encountered in almost all simulation 
studies- this time step value was found to be not too small. The system was 
thermostatted at the ends of the MD cell only, but very similar unpublished 
results with less variable fluctuation was obtained by thermostatting each layer 
with strict conservation of momentum during the thermostatting process 
which involves solving coupled equations for each thermostatted layer. But it 
was desired to mimic the actual experimental situation where the reservoirs 
occurred at the boundary location of the system. It is surmised that reservoir 
boundary conditions could play a pivotal role in determining the product out- 
comes of reactions sensitive to energy fluctuations. The term sole in the figure 
captions refer to the standard error. 



3.1 Equilibrium constants 

Two independent methods ( (i) and (ii) ) were used to confirm the thermody- 
namical results. 

(i) Time independent distribution sampling method 

In order to find the thermodynamic equilibrium constant, K eqi the following 
procedure was adopted. The concentration ratio, K c defined as 

K c = ^ (7) 

was determined as a function of average system density, p where the x's repre- 
sent number density concentrations of the species indicated by subscript for the 
reaction J5J. At very small densities, the system becomes an 'ideal' mixture, 
but as illustrated previously £Q , the limit of the potentials cannot be the same 
as the isolated potentials used in the MD calculations, since if this were the 
case, all the molecules would break up, yielding a net zero value for the equi- 
librium constant at the limit of zero density. As another project, it would be of 
interest to determine the limit at which the equilibrium regime breaks down in 
this thermostatted system, but there may well be technical difficulties involved 
in computations of very low density systems. The plot of K c = K c (p) is shown 
in Fig. ©. 

It will be noted that the errors in the K c ratio at low densities are of order 
20 times than that at higher densities, where relatively very accurate sampling 
is possible. What is observed at very low density is a '"saturation"' or limit 
effect for both rates (leading to fluctuations in K x in method (ii) and depicted 
in Fig. © and for K c for this method (i) given in Fig. ©. This allows us 
to determine K eq from taking the average value over a range of low p in the 
saturated range of p (here from p = 0.03 — .09) over about 13 values where 
theoretically K eq = limM — > oo J2i=i Kc,i/M over any saturated range. The 
resulting constant is 

K eq = lim K c = 0.061 ± .002LJ units (8) 



p^0 



Knowing this value, we calculate the activity coefficient ratio, $ e = for the 
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Figure 2: Variation of concentration ratio K c with p, the system number density 
at LJ temperature T* = 8.0 with sde — 3 at p = .03 and sde — 50 at p = .9 
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Figure 3: Low density fluctuation of concentration ratio with system density p. 
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other densities using 

K eq = K c -^- = K c <S> e . (9) 
The ratio of activity coefficients <& e is shown as a function of density in Fig. 10]). 
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Figure 4: Variation of $ e and (&k with /?, which is the system number density 
at LJ temperature T* = 8.0 where there is empirically complete coincidence 
within experimental uncertainty. 

It is evident from the plots that the mixture is highly non-ideal, which may be 
expected due to the large differences in the LJ energy well for the molecule and 
the atom (see Fig.Q. 

The separate activity coefficients may perhaps be derived from cycle studies 
of equilibrium states that include the reference state at infinite dilution. The 
derivation might require a series of very elaborate and detailed computations. 
However, to date no clear theory nor example has been provided as to how this 
might be achieved. 

(ii) Kinetic sampling method 

From the way the algorithm was constructed for molecular formation, the 
molecularity of the elementary reaction is 2 leading to a single second-order 
reaction of formation, and for the dissociation of A2 , a first-order reaction results 
since the molecule can only exchange kinetic energy with all other particles 
within the system without further reactions to the dissociation limit. There 
are many definitions and conventions for describing the overall rate of reaction 
and the individual rate constant. Here, the forward rate constant and the 
dissociation rate constant k°_ 1 are defined as the value of the respective rate 
constants as p — > where the unsuperscripted fc's refer to the rate constants for 
non-zero finite p. The overall rate of reaction r may be written in terms of the 
experimentally determined forward rate r% defined as r\ = k\x 2 A and backward 
rate r_i defined as r_i = fc_iXA 2 where r = n — r_i = k\x\ — k^iXA 2 - In this 
definition of rate, no normalization by the stoichiometric factors of the chemical 
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reaction are used. At equilibrium r = 0, and so 



^ = -£L. (10) 

X K— i 

The ratio of rate coefficients is the concentration ratio which is this time written 
as K x and is expressed only as a kinetic ratio where 

(id 

To verify the above we plot 

ri/x\ = Q = ki 
and 

r-x/xA 2 = R=k-i (12) 

against the density p and extrapolate to zero density to determine the equilib- 
rium constant. The rates were calculated independently from the program by 
monitoring the number of bonds formed or broken for each time step 6t* and 
averaging this quantity over the 10M time steps. The plots of Q and R at low 
densities (where saturation is observed) are given in Fig.© As with K c , a sat- 




Figure 5: Very low density fluctuation of R and Q. These values were averaged 
to yield the low density limits. 

uration effect is observed amidst the evident fluctuations at very low densities 
and the zero density (superscripted 0) limit for Q and R is defined thus: 

lim T M O 

Q° = fe? = M - oo 



M 

lim V M R 
= lvi — > OO — 

1 M 



( 13 ) 

R° = k°i = M -»• oo • 
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The subscripted Q's and R's in l|13|) are the values determined in the saturation 
region interval of p where 0.03 < p < 1.1. In this work M = 13. If the simulation 
is to be consistent and at equilibrium, then K c = K x = Q/R which is precisely 
observed in Fig. JSJ. The individual variation of R and Q are given in Fig. (Q). 
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system p 



Figure 6: Variation of K x with system density, with the uncertainty of sde = 50 
at p = 0.4. 

The results for the low density limits are as follows: 

lim Q = k° 1 = 0.870 ± 0.006L.J. units (14) 

p^0 

lim R = k°_ 1 = 14.32 ± 0.1L.J. units (15) 

and their ratio is 

ki 

K° c = lim — - = 0.0610 ± .006L.J. units (16) 
p^o A;_i 

An excellent agreement with the equilibrium K eq is found, i.e. K eq = K® = K®, 
where the method used for the determination of the equilibrium constant differs 
from the static distributive sampling of method (i). 

3.2 Survey of some standard form for elementary rate con- 
stant form 

An elementary reaction has been defined irreducibly [1 11 p.l] as due to the 
actual scattering of particles where "'..Chemical reactions occur by the collision 
of molecules, and such an event is called an elementary reaction for specified 
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Figure 7: General variation of Q = k\ and R = k-i with p, the system density 
at LJ temperature T* = 8.0 with sde(Q) = 50 and sde(R) — 3 at the indicated 
points. These values are at typical system densities which were not used for the 
extrapolation process. 
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reactant and product molecules" ' where if there are two or more elementary 
steps involved, as in the Michaelis-Menten reaction, then a " 'pseudo-elementary 
reaction" ' is implicated ^2 p-3].The rate coefficient (lit p. 7, Chap. 3] is written 
k(T), "' and is generally a function of temperature, and frequently of T only."' 
The onset of perturbations ^2 eg- Chap. 11, Oscillatory Reactions ] are due to 
species effects not connected to k(T) where the rate v of the elementary reaction 
is written v = A:(T)[^4] Q [_B]^... and the entire field of reaction mechanisms does 
not question this representation. The form was established over a century and 
a half ago and follows from the law of mass action Chap. 3] which was 
suggested by persons which included C.L. Berthollet (1799), M. Berthelot and 
P. St. Gilles (1862-63), L. Wilhelmy (1850), H. Rose (1842), A.W. Williamson 
(1859) and C. Guldberg and P. Waage (1864-67) where they all in different ways 
generally inferred that at equilibrium, 

k + R 1 R 2 .... = k_P 1 P 2 ... (17) 

where the i?'s and the P's were effectively the mass concentrations of reactants 
and products and where the fc's were " 'coefficients of velocity" ' which were in- 
dependent of concentration. This form has been maintained to the present times 
no matter how complex the reactions (e.g. possibly involving multi-step clus- 
ter and activated thermal electron transfer J3| or complex chemical oscillations 
|14j . In the preface of the text devoted to non- linear chemical kinetics ^J, ^ 
was surmised "'..only first order processes escape the non-linear net, and even 
these get caught if there is the slightest departure from isothermal operation." ' 
On the other hand, in the dimer reaction described here, the first order A 2 — > 2A 
process is shown to be non-linear in the most ordinary circumstances. Complex 
concentration effects on the rate constant have been structured by postulating 
the concatenation of several elementary reactions, meaning that these reactions 
are pseudo-elementary according to Ref. ^I]> where each elementary reaction 
rate constant is strictly concentration independent in all the references hitherto 
encountered, and which therefore can be taken to be a universal definition. In 
Eyring's transition state theory, a scheme 

Reactants f± Activated Complex — > Products (18) 

is postulated for the pseudo-elementary reaction l|18|) above ^1 eqn. 4.3, p. 125], 
leading to a rate R 

# = fc'(T)^...CiC 2 ... (19) 

7 T 

where fc'(T) is strictly temperature dependent, and the 7's are the activity coeffi- 
cients; k'(T) is derived from a product function involving equilibirum constants, 
and strictly elementary rate constants, implying no other variable dependency 
other than T. Likewise, pressure dependency of unimolecular and association 
reactions 1995 edition, p. 121-142] are explained by Lindemann-type mech- 
anisms involving complex equilibria, which introduces a pressure dependence 
on the pseudo-rate constant ^3 1995 edition, p. 138, Sec. 5.11, Association 
reactions]. Ref. ^3 Sec. 5. 4, p. 186 ] give examples of composite elementary 
reactions, such as ^3 eqn. 5.90, p. 191] 

A + M % A* + M, A* ^ A^ -> P (20) 

fc-2 



11 



where various quantum statistical approximations are applied to each elemen- 
tary step involving the individual kt(T) rate constants to derive the overall rate 
constant which is first order at low reactant concentrations and second order 
at high concentrations. As stated elsewhere, Kosloff has gone beyond the TST 
pseudo-elementary reaction theories by developing a straight trajectory model 
|171 p. 187] which superimposes trajectories with energy grids. More recent 
methods, such as developed by Miller [SJinvolves going beyond TST where an 
exact trajectory is calculated " 'which is no longer a transition state theory" ' 5 
p. 387]. However, in these exact treatments, the isolated participants only are 
included, e.g. p.eqn. 46, p. 402] in 

H + 2 — >OH + (21) 

where the CRP (cumulative reaction probability) is calculated for the 4 atom 
system (6 — D calculations). The microcanonical rate constant k(T) for a bi- 
molecular reaction A + BC — > AB + C is defined such that k(T) is the canonical 
ensemble rate coefficient in the expression — = fc(T)L4][i?C] and k(T) is 
determined from an integral involving the CRP. A direct evaluation can also be 
made without CRP's p. 408, eqn.50] where the direct calculation involves a 
flux operator T which is related to the system Hamiltonian and trajectory which 
in all these treatments does not include the rest of the environment variables. 
A recent example of such direct calculations in given in JHj for gas phase reac- 
tions, where the flux operator is a type of commutator T = j?[H,9(s)], where 
H is the system Hamiltonian and s are the trajectory coordinates with 9 being 
the heavyside function. Reference jf]] is an excellent resume of the most up to 
date prominent methods available, all conforming to the k = k(T) assumption 
for elementary reactions where no direct relationship with the activity of the 
species has been deduced, as will be attempted here. A brief survey of the re- 
cent past also does not yield any exceptions to the standard definition p. 36, 
p.109, eq. 3.25, p.239, eq. 7.38]. Espenson (23 p.l57]too, writes k = k{T) for 
elementary reactions. In developing the standard Bronsted-Bjerrum equation 
for the reaction A + B — > P (products), he breaks down this reaction into a sum 
of elementary reactions [201 eq.9-22] 

A + B ^ [AB]t -> P, (22) 

where he infers [201 P-204], in keeping strictly with the rate constant defi- 
nition that "'...The rate is taken to be proportional to the concentration (not 
the activity !) of the transition state"', thereby deducing the standard form [201 
p. 204, eq.9-27] for the Bronsted-Bjerrum equations as 

k(T) = k ref ^- (23) 

where in the present simulation, the reference state is the vacuum. Other trea- 
tises do not differ in interpretation when they write |21l p. 3]" '...Typically, rate 
constants are independent of concentration but they may be quite sensitive func- 
tions of the temperature." ' This independence of concentration for elementary 
reactions is maintained elsewhere, e.g. [221 p.l3,eq.Tll],|23l p.2]and |241 p. 35, 
eq.2.39]. 
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3.3 Reactivity, activity and availability coefficients and ra- 
tios 

Write 

R = k°_ l7A2 = 

where j' A and 7^ are denned as '"reactivity"' coefficients. In this system, the 
substrate in which the chemical reaction occurs is the vacuum state with a 
defined zero of energy relative to the component species. Obviously, j A — 

and 7 yl2 = j^- are computable since all the other terms are. These coefficients 
are graphed in Fig. JHJ. We make the following observation: 




Figure 8: Variation of reactivity coefficients j A and 7^ with p, the system 
number density at temperature T* = 8.0. 

Observation 1 j A and J A are not in general unity, especially at higher p 
values. 

According to the prevailing theories over the centuries, the rate constant proper 
is independent of concentration for strictly elementary reactions. Activity coef- 
ficients are only introduced via an equilibrium constant as a result of composite 
schemes where some pre-equilibrium is postulated, the chief example being ac- 
tivated complex theory, and the earlier and related ionic reactions based on 
the Bronsted-Bjerrum equation. In the present scheme, pre-equilibria has been 
ruled out, because the non-reversible, cyclical pathway of the reaction does not 
allow for such a situation. The results here show that the conventional definition 
of the rate constant for elementary reactive processes is incomplete, unless the 
7 reactivity factor is included, as summarized below. 



13 



Lemma 1 Given that the temperature- only dependent portion of the rate con- 
stant k B is always present, then its product with the reactivity coefficients are a 
necessary and sufficient condition for the rate constant to be complete. 

Proof. Denote by [X] the presence of factor X and V A the logical 'or' and 
'and' conjunctions respectively with -> the negation. Let the rate Sb for species 
B be given by 

S B = k° BlB (25) 

where j B are the reactivity factors (they may be product terms as for our A 
species). Then [7s]A[fc^] Sb is complete since form l|25|) leads experimentally 
to an exact description for all p values. The converse is 

- ([7b] A K]) = A1b\ V Ak B \ => ^B. (26) 

Since by hypothesis [k B ] is present, then [^fc^] is false and so if r y B is not present 
implies Sb is incomplete (which is impossible) • 
Fig. ©provides for the following: 

Observation 2 $fc = -^f- is coincident with <I> e within computational accuracy 
for all physically feasible p. 

This is an remarkable and unexpected result from which further properties 
issue forth such as: 

Lemma 2 The general solution to Observation^) is 

l\c{p) = (27) 
IaAp) = 7a 2 - 

where ja and ja 2 are ^ ae activity coefficients. 

Proof.The reactivity coefficients for any species X is known i.e. "f x = "fx(p) 
and r )A 1 lA 2 exists by hypothesis (Observation ©. Then 3 di such that "f\di — 

•yj and similarly r )A 2 d i — j Az .Observation© implies = l^j for all p or 

d i = di for every p,. Define the function c : pi — > d, or c = c(p) • 

If c(p) is a more complex function than unity, we may rescale our apparent 
activity coefficients such that we may relate the kinetically derived 7 reactivity 
coefficients with the 'actual' 7 coefficients determined from equilibrium distri- 
bution studies, so that 7^ = 'jac(p) 1 ^ 2 ,"/ A ~ lA 2 c(p) are used instead of the 
7 activity coefficient in the chemical potentials for the system. That is, if the 
actual potentials are scaled, we wish to determine whether j A and j A may be 
used as equivalent activity coefficients when ja and ja 2 are scaled by functions 
c(p) 1 / 2 and c(p) respectively. As a matter of interest, in the sections that follow, 
it will be shown that the chemical potentials refer to purely isothermal reversible 
work terms. 

Lemma 3 The scaled activity coefficients yield equivalent expressions for the 
Gibbs equilibrium criterion for the same Gibbs standard chemical potentials, but 
the Gibbs-Duhem condition demands c(p) = constant. 
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Proof. Writing the rescaled chemical potential for any species p as /i+ = 
fx* + i?T In Xp7 p and /j. p = /i* + RT In x p "f p for the unsealed potentials, then 
application of Gibbs' criterion for the scaled variables yields: 2\i\ = => 

2/ia = f-A 2 by algebraic cancellation and 2fiA = /M 2 = ph\ 2 by algebraic 

addition of lnc on both sides of the equation. Hence either set (/x or n + ) of 
potentials may be used to determine the equilibrium point where 7 p represents 
the activity coefficient for the /i p chemical potential set. The Gibbs-Duhem 
equation for the scaled system at constant (P, V) is 

J2^dnt=Q (28) 

i 

where 



d/i+ = dfi* p + RTd In x v + RTd In lp + RTd In f p (c) . (29) 
Comparing i|28[) with J^. Xidfii = for the set, we derive 

x A RTd\nc 1/2 + x A2 RTd\nc = Q (30) 
or 

—^-dhxc = XA 2 dhxc (31) 

which is generally a contradiction unless c(p)=constant, which implies c = 1 
since lim p ^o c{p) = 1 • 

Corollary 1 From Lemma pfy. either set (\i or^i + ) may be used to determined 
the equilibrium point, but there is only one specified activity coefficient set {"f P }- 

Corollary 2 If it can be proved that for any rescaled activity coefficient^, there 
can exist only one equilibrium concentration ratio satisfying the Gibbs equilib- 
rium criterion and the Gibbs-Duhem equation for the same chemical potential 
functions, then c{p) — 1 and = jj . 

Corollary 3 The activity coefficients are unique in that if they are rescaled, 
then c(p) = 1 only. 

Corollary 4 From the invariance of the equilibrium constant due to the Gibbs 
criterion, the rate constants between two states 1 and 2 due to catalytic activity 
are related as follows: 



kl \ fk-i\ (~fM\ =1 (32) 



k-lj 1 V h J 2 V T A ) i V7A 

One may not expect the 7's to vary between two equilibrium catalytic states 
by mutual cancellation but ea . I)32|) is the general expression. Given the same 
mechanism, changes of rate can only be affected by changes in the isolated 
intermolecular potentials, or the net potentials due to the sum of all relevant 
interactions. Both these cases seem to imply minute changes to the species type. 
Provided that these changes do not affect to first order AG° for the standard free 
energy change, ea. (|32|l should obtain approximately. Since Q, R, and the 7's are 
functions of p, we can write equations i|12|) in the form 
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R = fc ^(X>7i^7A 2 (33) 

If R and Q are independent variables, it is not immediately obvious why the com- 
mon factor c(p) — X)i=i a il\ = Si=i ^7a 2 mus t have the same non-constant 
value at each p if we should vary the temperature: hence one might expect that 
c(p) = 1 : if this is not generally the case, then there is embedded within the MD 
potentials a yet to be clarified link between the forward and backward rates. 

The plot of the reactivity coefficient is given in Fig. JHJas a function of aver- 
age system density. Both coefficients extrapolate to unity. While the reactivity 
coefficient of A is larger than unity, that of A2 is smaller than unity, and one 
might expect such behavior to arise from the differing potentials and net surface 
area to volume ratio that exists for these two species. The reactivity of species 
p,a p is defined by a p = Jp-Pp where p p is the concentration (e.g. number den- 
sity) of species p. A plot of the variation of these quantities appears in Fig. JHJ), 
where the densities refer to the actual number density of the concerned species 
and not to the general system density p. The looping of the curve for A is very 
interesting because with increasing system density, the equilibrium constant ad- 
justs itself to accommodate an increasing reactivity coefficient for A, implying 
at constant temperature a species density reduction, which explains the looping 
back behavior of the reactivity for A. At the current state of development, one 
must use physical arguments to determine whether the reactivity and activity 
coefficients are identical or not. It is surmised that it is not from considerations 
of the estimates of the LJ fluid activity coefficients relative to the boundary 
conditions imposed. It is estimated that both the dimer and atomic activity 
coefficients either increase or decrease with system density, and therefore can- 
not equal the reactivity coefficients where the trend for dimer and atom are in 
contrary directions. It is then surmised that the activities both increase with 
system concentration by comparison with the LJ fluid for the atomic constituent. 
Before examining the thermodynamics of the single component monomer A, we 
digress to an outline of a theory of internal species equilibrium that is of use to 
describe species states along a reaction pathway. 

3.3.1 Internal equilibrium of product species 

It was conjectured that a subsystem described by non-canonical coordinates 
would have a Boltzmannizcd energy distribution for the kinetic energy if the 
energy interaction of the coordinate were due solely to external forces (due to 
the reservoir). This conjecture is extended to the external potential here; in 
fact the Debye-Huckel (DH) and extended theories of electrolyte solutions all 
make this assumption |25| . Under this assumption, a particular species Z (being 
either product P or reactant R species in the equilibrium R P) can exist in p 
varying energy states Z±, Z2, ■ ■ ■ Z p as enumerated by an appropriate algorithm. 
The enumeration may be with respect to energy states in terms of a coarse 
grid of magnitude AE; such techniques have been developed in recent times 
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Figure 9: Variation of reactivity a a and aA 2 with pa and pa 2 , the actual number 
density of the species concerned at LJ temperature T* — 8.0 with sde — 10, 40 
respectively for the uncertainty in a a and dA 2 . 



by Kosloff [171 p.l87]who superimposes reaction dynamics trajectories with a 
gridded "'energy range of molecular encounter AE ...typically in the range of 
— 10 eV."' And It will be further assumed that a convenient algorithm ex- 
ists for this purpose. The example given here is for reasons of illustration, to 
show that the activity coefficient may be determined at any point if a partial 
enumeration were known, and to extend the DH theoretical assumption. The 
science of enumeration has made much progress in recent years, especially in 
the characterization of species types in mathematical chemistry, whose pioneers 
include Balaban 26 . The techniques developed for chemical structures could 
conceivably be extended to these cluster types. The temperature of each of the 
species in the enumeration is approximated as the same as the system temper- 
ature T determined from = / Ji_ ^ where it has been shown that they may 

differ slightly since they are not the canonical coordinates Hence this theory 
is an approximation. Then the chemical potential of each of these states 
share the same standard state p? z where \iz i = + kT In Ciji with Cj and ji 
representing the concentration and activity coefficient respectively and for the 
Z species as a whole, its concentration ctot will be Ctot — J2i=i c i an d the bulk 
activity coefficient will be denoted as jtot- The species equilibrium 

Z 1 ^ Z 2 ^ Z 3 . . . ^ Z p (34) 

implies a form 

kTlriCgmjgm = kT In c l ^ l = kTlncjjj (j ^ i) (35) 
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where c gm and 7 9TO will be determined. From ((35(1 . let atj — Ciji, then a,i — 
otj = a' since p,° z .(T) = p.%{T) , and further define 



From (EH, 
leading to 

Define 



kT In | Y[ c i7i | = feTln A (36) 
kTplna' = kTln(cic 2 . . .c p ) + kT In (71 72 ■ • - 7 P ) 
fcTW = fcrin{ Cl c 2 . . .c^} 1 ^ + fc7Tn{ 7l 7 2 . . . 7p } 1/p (37) 



lnc gm = ln{cic 2 c 3 . . . c p } 1/p 

In 7am = M717273 ■ ■ -c p } 1/p . (38) 
Then c 9m = {cic 2 . . . c p } 1/p and -y gm = {7175. . . . J p } 1/p . Further Ogives 

kThiCg m 'yg m = kThia' = kThijiCi. (39) 
It will be noted that © 

can refer to a partial (incomplete) enumeration where 
the i enumeration need not be in the set if the Gibbs equilibrium criterion 
obtains for the i th species in equilibrium with the other species. 
What is the relationship between the different c and 7 coefficients? 
For a complete enumeration over an energy grid, each of which has an energy 
difference of magnitude AE, the average potential experienced by the species X 
is W, so that w = kTlnjtot by definition, where [iz = n% + kThictotltot- The 
probability of state i, Pi is pi = ct/ctot and under the the canonical ensemble 
assumption above, pi — cx p ~ e '^ wnere z = X)f=i ex P — e ifi is the partition 
function of the distribution with (3 = l/(fcT). By taking averages, the above 
yields 



or 



fcTln 7iot = kT^pilnji (40) 

2 = 1 

P 

ln7t oi = ) — Inji (41) 



i=l 



Hence, (|4*T1) leads to 



ctot In Jtot = y^Cjln 7i (42) 



Defining a departure from bulk value 7' where 7, = jtotJi, then 142(1 can be 
written in the form 

by expansion of terms. 

If the degenerate terms are all collected together, the canonical distribution 
can be written pi — 9i cx P^l 3ei , (Other modifications are straightforward and 
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here Z = X]f=i 9i ex P ~P e i)- We assume p non-degenerate energy states. Since 
ji = exp (3ei , then the probability term is pi — -^g and so 

In = In g l - In Z - In (44) 

so that (|4*T)> becomes 

In 7, 



p p 
hx-ftot = TTz ln71 = ~z ^ 9i 



and for gi = \ there results 



(45) 



ln7 fot = ^lnn( 7 - /7 -). (46) 



What can be said of jtot in terms of c gm ^ gm = a where <|39[) gives c gm ^ gr 
a = 7iQ? From these definitions and (|42(1 . there results 



^^^^^^^/^ 1/7 A (4?) 
Ctot i V ' 

2—1 

and comparing (|47|l with (|46|l implies 



Ctot 



For degenerate systems, let ai = gi/ji, then 



(48) 



(49) 



1 v 

Then 147|l and l|49|l leads to a general coupled equation 

P P / cgmTgmA 

ik i/z =in ct ° tTi J - (50) 

i=l 

Forcing independence of the 7j is equivalent to equating exponents, leading to 
the result 

2 C tot 

Similarly, (|42[1 and the definition of a leads to 

p 

ln7tot = In a - J^Pi lncj. (52) 

i=l 

Using the definition for Cj in terms of pi, we derive 

a P 

In 7tot = In ^ Pi In (53) 



Cto* . , 

2 = 1 



which has an interesting entropy-like contribution in the probabilities. Finally, 
the maximum available potential energy for any one species 7$ must obey ele- 
mentary relations such as 7 9m c gm = 7^0^ and ^j 2 - = 1, implying that one can 
determine the 7 terms if the concentration terms were known at any volume 
element of the chemial trajectory. 
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3.4 Thermodynamics of fluid and discussion of the c(p) 
scaling function 

LemmaQ implies that since the activity coefficients are unique relative to scal- 
ing, an attempt to estimate c(p) can come only from physical/kinetic experi- 
mental considerations, where this function should not contradict ancillary data 
and estimates. One can expect the activity coefficient of the A atom to be 
approximately that of the LJ fluid at the same density and temperature at low 
dimer concentrations. A rationale must be provided on how to estimate this 
quantity from available data. Some data exists |57| |2H| for the L J fluid up to 
T* = 6.0. At such temperatures, the fitting is highly sensitive to the 7 coeffi- 
cients, unlike what was reported 1281 sec. 4, p. 607] (presumably for much lower 
T values, where 7 was regressed for 1 < 7 < 7 with low correlation of the mini- 
mum). The parameters for the various estimates for the monomer activity was 
scanned from the x(i) values in Table 4. of [22] and fitted into the equations 
found in [2Hj- The residual or excess Helmholtz free energy in reduced units, 
which is the difference between the ideal gas value and the actual value is given 
by [2D eq. 5] 

^■ = E^ + E^ (54) 

i=l i=l 

with the pressure having the expression 

P* = p*T* + E a lP <^ b lP * (2i+1) (55) 

i=i i=i 

with F = exp(— 7P* 2 ), 7 being a non-linear adjustable parameter. The excess 
free energy G* is given by 

G* — A* + — T* . (56) 

P 

The reduced units variables are denoted by * superscripts, and the and hi 
coefficients are given as complex polynomials in x(i) and T* in the tables in 
|28| . It was found that at T* = 6 , a good value for 7 (to the nearest digit) 
was 2.7 which yielded the fitting pressure of f2.7931, close to the experimental 
value of 12.43(2) in the phase diagram. The simulations here were off-scale at 
T* = 8.0, but nevertheless the value of 7 given above was adopted for estimating 
the excess free energies at the simulation temperature as well as the atomic "/a 
coefficients according to the theoretical model below. 

A circular reasoning has been detected in the definition of multicomponent 
Gibbsian free energies |29l p. 706, Sec. 5] and the associated chemical potential, 
where it was pointed out that the TdS term defined by Denbigh 30 sometimes 
referred to a heat input term for multicomponent closed systems, and some- 
times not, leading to an undecided paradox, but where the chemical potential is 
a pure work term; in the new exact thermodynamics [291 p. 707, after eqn(30)], 
the chemical potential is a mixed heat-work term. In this work, we shall con- 
form to conventional descriptions, thereby reinterpreting the physical meaning 
of some of the standard expressions. It will be argued here that relative to 
the concentration where the chemical potential is written /i^ = /i° + kTliiCi^fi, 
and where this term represents isothermal work done on the system, the excess 
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Helmholtz energy A* seems to refer to the 7^ activity coefficient (here analogous 
to the '"fugacity"') coefficient and not the excess Gibbs free energy G* for the 
reasons that follow. For a single phase at constant temperature, dG — Vdp 
and for a perfect classical fluid AG = / Vdp = — J pdV where the work loss 
on the system is equivalent to AG but this is not true of imperfect systems. It 
will be proved from the Kelvin-Clausius statement that the external work (if it 
were the sole work source) takes into account the internal intermolecular forces. 

The chemical potential tii = ( M^- ) , despite being defined as EH P-79] 

"' the amount by which the capacity of the phase for doing work (other than 

work of expansion) is increased per unit amount of substance added...." ' is often 
written in the form 

fU = £{T) + kT\n Clll (57) 

for substance i. The definition by Denbigh may perhaps be reflected in the view 
of some others who suppose that even for single phase systems, it is believed that 
G/N = fii 31, pg 83]. Indeed, for a single phase system, some have identified 
the chemical potential with the Gibbs function as U — TS + pV = /in = G with 
n being the amount parameter |31l Sec. 7. 3, p. 83]. 

This would suggest that for multicomponent systems, the free energy differ- 
ential written 

dG = -SdT + VdP + Vidrii (58) 

i 

implies that the chemical potential refers to a contribution to the free energy 
to the state at constant p and T |S2 p. 126] with the material addition drii. 
Then from G = ^ n iA*i |32l p.255,eqn.9.23] , the one component form suggests 
that = fii or d^'"'- 1 = ^ or dCfcTinc^) = for all species i which is 
undemonstrated. Furthermore, for imperfect fluids J VdP ^ — J PdV generally 
and so does not appear to constitute a pure external work term for single phase, 
single component systems; if it were nevertheless the case, then 

G eXt i = Greaij — G ic i,i = fJ?(T) + kT In a-fi - fJ?(T) - kT In c, 

^7, = exp(^). (59) 

To estimate the activity coefficient of the LJ fluid, one must refer to the 
excess work done on the particle 5w ex ^, where conventionally, the activity co- 
efficient 7, refers to this work. [251 Chap.l,p. 4-8] where kT\n-fi — Sw ex ^; 
it will be shown that the term kT In Cj refers to the external work done on the 
system to the state concerned, and this work is also equivalent to introducing 
unit amount of substance i from the standard state condition of unit activity 
coefficient. The work done for a perfect fluid on the system is 

Sw = -NkTln— = NkT\n— (60) 
Vi ci 

between states 2 and 1. Under the typical convention c\ — > 0, 71 — + 1, the 
standard state c\ leads to a singularity in this limit. Write SciMi — 1 for any 
small value of c\ of amount 8c\ when 7, — * 1. So, Mi — > 00 but for any small 
8ci, it is a large finite number. The form of the chemical potential, if rescaled 
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to conform to these limits must be of the form 

rf-tfrc + H-b.^. (61) 

But whatever the scaling, the form and value for these variables match, i.e. 
/i; = /4 where //? refers to a standard chemical potential of isolated substance 
i as Ci — > since the physical meaning of /ij must be preserved. Thus, 

/z°(T) + fcTln ^ = M °(T) + fcTln M« + fcTln (62) 

so that 

V° i '(T)= t i i (T)+T.Q i (63) 

where 

Q, = k In M». (64) 

Let substance i, also termed substance A be the basic constituent from which 
all other substances are formed. Other substances are built upon this primary 
substance, where in general 

A n + A m ► A m + n (65) 

and the subscript m refers to the number of units of A making up substance A m . 
For instance for the reaction nA ^± A n , the standard state chemical potential 
^a„ CO f° r A n is defined as 

n\ n {T) = n^ A + 8^ An . (66) 

Here, Sfi A refers to the work required to form A n from n units of A. For all 
reacting species A q , a common lower end concentration limit is imposed where 

5c An = Set = Sc A = Sc Aq (67) 

for all n and q. Because of this, we can write 

^ A ' n {T) = n. l i A l + 8A n (68) 

where 

f i A , =»° i (T)+TQ i (69) 

with Qi = kin Mi. But from 

nv% = nfj#(X) + nTQi (70) 
and from (|68|l the following results 



MX, = ^(T)+nTQ l + M„ (71) 

= n/i°(T) + nTQ An + SA n (72) 

But l|66(l yields by comparison to H72J) 

4=^ + ^. (73) 
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The equilibrium criterion gives 



M M 

-AG = - AG ' = - /J} Vj = - E = kT m K e (74) 

3=1 3=1 



and ijZBt gives 



A I M M 



Y flv, = E fyi + E ( n i T( 3) ^ (75) 



Now , Tij is the number of elementary species i that constitutes j, and this 
quantity is always conserved (from a chemical point of view) and T and Q are 
common, leading to TQJ^j n j v j — 0- This verifies 1|74[) and infer that the 
same equilibrium point is reached. It follows that we may ignore the terms 
that cancel off if only equilibrium problems are of significance, and write the 
chemical potential for equilibrium problems as 

^(T) = M ?(T) + fcTln^ (76) 

where the potential here is a measure of the work done on a particle i in isolation 
where the singularity of the point at oo (c$ = 0)has been removed. The above 
establishes the fact that the form n"(T) = /J>°(T) + kT In aji is an isothermal 
work term with removed singularities at zero concentration, consisting of the 
work used to overcome the internal forces and the external ideal work. 



3.4.1 On Helmholtz heat- work interchange 

An approximate value for the multicomponent activity coefficient may be de- 
rived from the Helmholtz free energy for single systems by considering the work- 
heat transformations. The single component Helmholtz free energy A is given 
in standard form as A — U — TS, where its differential dA = dw — SdT implies 
that at constant T , the work done on the system dw is the so-called external 
work (such as the P — V work connected with compression) where for fluids, 
dA = —PdV = dw. For fluids, 3 basic categories ((i)-(iii))of compression from 
state at concentration 8c\ to c (where the lower concentration is at Sc±) must 
be examined as follows: 

(i) ideal fluid compression. 

The reversible work of compression is Sw = kT In and let the heat ab- 
sorbed be — Sq. The energy change between the states (subscripts denote the 
state for the variable concerned) AEi = E c — Eg Cl = Sw — Sq. 

(ii) compression of gas by by manipulation of internal forces. Here two separate 
sub-parts are considered where (a) is generally applicable and is in the standard 
format of a conventional thermodynamical system and (b) develops an inter- 
nal potential description with control of the particle potential. For both cases, 
divide the path into interval segments with states S indicated by subscripts : 
So (at Sc±), Si, S2 ■ ■ .S„(atc). For (b),we suppose that the external pressure P 
of the imperfect gas obeys P > Pid, due to the repulsive intermolecular forces 
relative to the perfect gas with pressure Pa for any concentration q for state 
Si. For a segment for which the opposite inequality holds, the same conclusion 
of the theorem obtains and so a general path may be broken up into subsec- 
tions with one or other of the inequalities obtaining (where it is assumed that 
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the pressure is continuous over the path), implying the general validity of the 
theorem. Method (a): Start with a perfect gas and compress reversibly to S n . 
At S n , we "'charge up"' the gas by introducing a potential with an interaction 
parameter e say. In Debye-Huckcl theory of solution activity, where the charge 
of the ions assumes the variable e; in a LJ type potential system, a possible e 
parameter for the potential /z(r) is /i(r) = e (-^ — rrr). For any potential and 
given particle coordinates, we can in principle compute the total potential of 
interaction. For each particular value of e we can calculate the mean energy of 
interaction per unit increment of e. Integrating this yields the average energy 
of potential interaction as WJ^. During this charging up process, heat energy 
will be absorbed by the system of amount —Sq' m during the manipulation where 
the system is at a fixed temperature T. For the entire compression process, the 
total heat absorbed is 

-Sq' = Sq' m - Sq (77) 
Then the total change in energy is 

AE(ii)=5w + wJ^-6q' (78) 

where Wf a is perhaps a novel diathermal work term involving internal potential 
variables; the standard methods hitherto used seem to suppose no heat transfer 
during the charging process and possible fundamental flaws in the thermody- 
namics might be introduced as a result. 

(iii)The third process is reversible isothermal compression of the imperfect gas 
by doing P — V work from So to S n at concentration c,where — 5q c of heat is 
absorbed by the system. The total change of energy will be given by 

AE(iii) = Sw c - 5q c . (79) 

Since the final and initial states are the same, then the First Law implies 
AE(ii) — AE(iii) but nothing thus far can be said about the Sq's, i.e. Sq' 
need not equal 5q c . Suppose in fact that the heats are not the same Sq' ^ Sq c . 
A cycle of state transitions can be carried out as follows: 

step(ii) step(iii) 

For this cycle, the First Law yields 

Sw + WJ2 — Sq' — Sw c + Sq c = (81) 

leading to 

{Sw + wj:- 5w c ) = (Sq' - Sq c ) ^ (82) 

or the effect in this cycle is the NET conversion of heat to work about an 
isothermal cycle, which is in contradiction to the Kelvin-Clausius Second Law 
postulate p-101]. Hence, Sq' = Sq c •. 

For the above, we have extended the scope of the traditional understanding of 
the system as given to manipulation of the internal potential of force interac- 
tions. This result may be stated thus: 

Theorem 1 The heat absorbed by a system during an isothermal transition is 
equivalent to that which is absorbed by an ideal system and the heat absorbed 
when an intermolecular potential is introduced within the system. 
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Corollary 5 The reversible work performed during a system transition is that 
due to the work for an equivalent ideal system and the potential energy change 
due to intermolecular forces and therefore the excess work done relative to the 
ideal fluid is due to the intermolecular potential energy change of the system. 

Proof. The first part follows from Sw c — Sw + Wf a from (|82[1 and theorem ^ 
and the second from Sw c — Sw = Wf a • 

Corollary 6 The change in the excess Helmholtz free energy SA ex is equal to 
the change in the intermolecular potential energy for the system. 

Proof. By definition, SA ex = Sw c — 6w and this quantity is given by corollary[S] 
by Wf a u 

Corollary 7 The activity coefficient for a single component system 7; relative 
to the standard state is given by ji = exp ( ) . 

Proof. Since Sw c = kT In , Sw — kT In , and 5A ex — Sw c — Sw then 
SA ex = kTln^fi* 

In the dimeric reaction 2 A ^ A2, for lower A2 concentrations (lower p val- 
ues), the limit ja, pure (pure phase) — ► ja, reaction, the activity coefficient for the 
reactive system A is exactly as for the single component fluid because the force 
fields acting on A are essentially the same as for the pure component, except 
there would be clumps of A atoms held together by intermolecular harmonic 
potential whenever A2 is present, but these internal forces do not affect the 
external forces acting on the A atom. This identification and limit is used to 
choose the general form of 7 vs. p curve for atom and dimer, since the two differ- 
ent assumptions that are used leads in one case to semi-quantitative agreement 
of 7^4 for the reaction and 7a, pure (pure phase), and it is therefore inferred that 
the system conforms better to one of the two assumptions made to derive esti- 
mates of the activity coefficients. Clearly, only one of these must obtain from 
the previous deduction (Corollary 0) concerning the uniqueness of the activity 
coefficient. We can infer that the results from Fig. p5|l is the more consistent 
one by referring to the monomer situation, for which some data exists, and 
by estimating the monomer activity from the data and to then infer that the 
form and trend provided for the monomer at lower concentrations should be 
approximately that for our pure atomic species as deduced from readily avail- 
able functions. These functions were estimated for values of temperature up to 
T* = 6.0 and therefore are not adequate for higher temperature estimations. 
Nevertheless, even for T* = 6.0, the values follow exactly the same trends as for 
the results here at T* — 8.0, including the activity estimates all much greater 
than unity. The 7 estimates are derived from ean. (|59(l and the expression from 
Corollary @. It is clear from Fig. p()|l that the activity estimates show values all 
greater than unity. The system p for this figure and all others refer to p system, 
all in reduced units where p system = PA + %Pa 2 i where pa is the number density 
for the monomer (Atom) and pa 2 refers to the dimer species number density in 
this homogeneous system at equilibrium. Fig. \ 10(1 is computed for the situation 
PAtom — P system- Since what is being depicted is the external potential work 
done in bringing the atom to the system, then the forces acting on the atom 
can arise from either the dimer or the atoms within the system, and since the 
potentials are of the same LJ form for each nucleus, then we would expect this 
density to be the most appropriate one. If we just consider what the activity 
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Figure 10: Computed monomer-only activity from residue Helmholtz (A)and 
Gibbs (G) functions 

might be for the actual density of the monomer, ignoring its interaction with 
the dimer, then then Fig. l|ll|) gives the results from the simulation results and 
the fitting functions for G res and A res where the system concentration p = pa- 
Finally if we supposed that whatever energy that is available to the monomer 
is lost whenever a dimer is formed from the activation energy then the available 
energy per atom is given by Fig.H12[) . What is meant in this case is that let the 
residual free energy per particle be F r} where F r = G res or F r = A res .The vac- 
uum activation energy is 17.652. So the net energy left per particle is surmised 
to be F res = f ^a-i™52N A2 17 652 p A2 /p and 7Fres = exp(F res /(fcT)*) 

which is a measure of the available energy after dimer formation. This expres- 
sion is plotted in the figure mentioned above where Na is the number of atoms 
for that general system concentration in the MD cell. These figures are plotted 
to determine which would be the best approximation to the first order theory 
of elementary reactions developed here. 

3.5 Considerations in Hamiltonian systems 

Recently, it has been proved that the Liouville theorem, and in particular 
the Liouville space, of cardinal importance to both classical and quantum- 
mechanical theories of dynamics and statistical mechanics does not obtain 
typically in standard applications; and therefore that non-Hamiltonian theories 
posited using Liouville space, such as developed by Tuckerman et al. are equally 
suspect |34l Tuckerman et al. references within]. A work which does not im- 
ply a steady state but a disintegrating system which preserves a certain form 
of entropy had been described earlier without recourse to Liouville space and 
operators, and which considers the Jacobian of a transformation '"Non- 
Hamiltonian" ' systems has been variously described, including their relation 
to thermostatting using a pseudo-Hamiltonian that forces constant conditions 
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Figure 11: Computed monomer-only activity from residue Hclmholtz (A) and 
Gibbs (G) functions 




Figure 12: Estimated activities from the residual energy left after dimer forma- 
tion using density distribution data from the simulations. 
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such as pressure or temperature |38| . It should be noted that many of the syn- 
thetic thermostats were designed with the belief in the Liouville theorem. As 
such, it would be expected that equilibrium properties might be computed if 
the phase-space covered follows the canonical energy distribution for the coordi- 
nates. However, if the cause of the canonical distribution is not due to a smooth 
Hamiltonian trajectory, but to external random impulses, then it can be inferred 
that the real time behavior of particles described by such pseudo- Hamiltonians 
over any short enough time interval St would not follow the prescribed pseudo- 
Hamiltonian trajectory. Such an inference implies that the dynamical behavior 
of the ensemble of particles would not correspond to the actual possible trajec- 
tories with random perturbing forces. It then means that molecular simulation 
with such pseudo-Hamiltonians would not yield realistic trajectories which could 
explain real-time, non-averaged behavior over very short time intervals where 
random impulsive forces were acting; in particular, biological system function 
based on real time changes of coordinates cannot be accurately simulated by 
such pseudo-Hamiltonians if their main mode of dynamics involves random per- 
turbations, such as from a reservoir. Another presupposition concerns the static 
nature of a Hamiltonian system [38, p. 496] and the supposed effects of '"walls'" 
and constraints in the Hamiltonian: '"...The notion in phase space of Hamilto- 
nian system is similar to that of an incompressible liquid; in time the volume of 
the '"liquid"' does not change. In contrast, a non-Hamiltonian system is com- 
pressible. This compressibility must be taken into account when considering 
the generalization of the Liouville equation to non-Hamiltonian systems'". On 
the other hand, from the work presented in |34|.it may be inferred that for any 
Hamiltonian form 



one can consider the " 'walls" ' and other physical constraints of the system as 
perturbations to TIe where V(r) does not refer to the wall interactions; for per- 
fectly reflecting walls, for instance, a micro-canonical distribution results, the 
canonical one obtains for the single system. The system trajectory is Hamil- 
tonian over the time intervals Tj = U — ti-\ for the set of time coordinates 
{tQ 1 t\ 1 t2....t m , ...ioo} of interference by random energy interchange by the wall 
or/and thermostat at time U; here the walls and energy interchange achieves at 
least two functions (i) it defines the boundary dC and (ii)acts as a thermostat in 
terms of energy interchange and the change of system trajectory not determined 
by the Hamiltonian. If the energy levels arc quasi-continuous - which certainly 
is for this dimer system- then TLe will be distributed with a canonical energy 
distribution P(E) oc exp— ^TLe according to standard statistical mechanical 
theories p. 1-3]. For the Hamiltonian system given here, eas. <|3l6[) . the K.E. 
profile of (i) the " 'Atom" ' is confined to that segment of the potential curve 
where there is no bonding, and (ii) the product A2 to the molecular potential 
well, with (iii) the '"transition state"' point located at two different points 77 
and r;,. The coordinates of the defined quantities (i-iii) are not in general the 
canonical coordinates of the entire system Hamiltonian. These coordinates rep- 
resent transient species and the energy terms for these coordinates may or may 
not be canonical [HI, and if they have a Canonical energy (CE) distribution- such 
as the K.E. of the centre-of-mass of the dimer- then its apparent temperature 
defined as ^Mk =< K.E. c/m , > need not correspond to T s , the system temper- 
ature. It was discovered 151 that the mean K.E. of the atom A and dimer about 




(83) 
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the CM. followed the C.E. distribution, with apparent temperatures slightly 
differing from T s . This result is not unknown, since within a reacting system, 
apparently different temperatures have been experimentally detected 0U] f° r 
the different species in a plasma. Gibbs and his equilibrium criterion suggests 
that all species at equilibrium must have the same temperature; we observe that 
experimentally this is not the case, where relative to the CM., the temperature 
corresponding to the K.E. is close but not equal, i.e. Ta ~ T s ,Ta 2 ~ T s but 
Ta 7^ Ta 2 - Currently, it is unclear whether Ta — Ta 2 — T s as p —* but this 
assumption was used for a new theory of energy interconversion £Q. For an 
elementary reaction species i — > products at very low pressure or concentra- 
tion, the rate constant k® (T) (p — > 0) , is derived assuming that the work done 
against the interparticle force due to reactant Ai is derived from the kinetic en- 
ergy of the particles at a particular geometrical locus determined by the impact 
parameter and geometry of the potentials and the activation energy; all these 
factors f2 are almost always considered invariant in Q.M. /classical reaction rate 
theories since they are strictly mechanical properties, so that k® — fc?(f2, Ti )0 ) 
where it is conventionally assumed that T s = T, lG ; for this model, the minimum 
impact parameter is at r/ = 0.85 with the activation energy at 17.5153. At 
other concentrations, the reaction pathway system would have to be perturbed; 
this situation is not normally considered even in advanced treatments [2]; the 
perturbation suggested here are the energy levels of activation; there would be 
other effects, of first and higher orders which perturbation theory can provide. 
The object of the current work is to provide an outline of a theory with esti- 
mates of the more important effects which can account for Observation^ and 
Lemma(j2J)- A second order effect concerns the change of mean kinetic energy 
of the molecules from the zero state 8Wi : K.E. (meaning the state as the density 
p — > 0) where 

6W ilK . B . = K.E.(CM,T iiP ) - K.E. (CM, T i>0 ) (84) 

but would not be considered further; obviously this perturbation would affect 
the number of activated atoms/molecules. The general perturbation is from a 
reaction at zero density (almost always the form given in standard theories) with 
rate constant = fc°(£l,Ti. ) to that at any other density which is the actual 
equilibrium state (a.e.s.) ki = ki(fl,Ti(T s )); i.e. a mapping from these two 
states is required. In the ideal state formulation (for example SCT), absolutely 
no other environmental interaction <p e of particle/complex reactants i with j is 
envisaged up to the collision distance a. In the perturbed state, on the other 
hand, for r# > a, <f) e {r x ) ^ 0,x = i,j for the coordinate r x of particle x. We 
note that for the vacuum state for which is defined, 4> e = for any particle 
participating in the reaction. The work done 5wi >p in transferring a particle i 
from the standard state of zero density (identified with the rate constant fc°) to 
a.e.s. is trivially 

Sw hp — kTlnji (85) 

with 7j being the activity coefficient. The concentration term corresponding to 
the ideal external work does not feature since it is a common term for both the 
ideal and non-ideal states. In (|85(l . 7 can be determined from MD by calculating 
the mean potential </>ifor any system p, where 7; = exp Standard MD gives 
the same quantity by a particle insertion method |411 p. 349, Appendix C] where 
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kTlnji and ^ res has the form 



f-^res 



= -fcTln 



< kT in >3/2 



1 



( 



(fcT m ) 3/2 exp 



( 



-u t 



) 



) 



(86) 



with Ut being the instantaneous potential of the test particle interacting with 
other particles. An obvious generalization of eq.ljHHJ involves writing a Hamil- 
tonian of M primary particles with a boundary dC not incorporated into the 
Hamiltonian in terms of a potential, with instantaneous energy Ei, which would 
fluctuate when it exchanges energy with the boundary which incorporates a heat 
reservoir in the form 



Although the Hamiltonian Ij87(l is invariant in its form, the system dC deter- 
mines the density of energy levels and the relative energy levels; by equipar- 
tition, the average kinetic energy is invariant for all enclosed boundaries, but 
the potential energy is not, since it is possible to set dC such that the inter- 
molecular distance of particles i and j is |rj 3 -| < C and I could therefore arrange 
a lower bound by suitable compression and choice of the boundary dC such 
that V[rx,r%,....ri....ryi) > Bqc- The Canonical distribution for the Hamilto- 
nian would obtain if the degeneracy is very great, according to the principles 
developed in the Feynman development |39| . Thus relative to the zero density 
state, where the maximum \rij\ — * oo and min \V{r\, r^, ....rj....rftf)| = 0, the 
mean potential energy changes with the boundary (and therefore system den- 
sity); further, the mean potential experienced by dimer and reactant atoms 
would also differ with density changes; hence relative to the zero density reac- 
tion scheme, the additional potential energy present is available to overcome 
the activation barrier of the reaction. The total available energy is 6wi^ p , and 
for the ground state reactants, this energy would lower the apparent activation 
energy, and will be used to parameterize the change in activation energy; it will 
be noted that there might exist the possibility that not all this energy might 
be utilized in extremely fast reactions where "'molecular chaos'" and the ad- 
justment of neighboring particles is not rapid enough for the full utilization of 
this potential. Further investigations, both theoretical and experimental would 
be required to elucidate these remarks. Hence a flexible approach is required 
to interpret the 7's which may in some cases be a fraction / in value to the 
thermodynamical 7; the results here indicate / rs 1. 

The proposed models for particle interactions yielding chemical reactions for 
the forward (F) and backward (B) steps are given in Figs. I | 1311 41 below. These 
models enable one to calculate the activity coefficients where one model con- 
tradicts experimental observation whereas the other does not, in addition to 
corroborating the deductions of the activity from thermodynamical functions. 

The inverse <£> convention yields the activity coefficients given in Fig. (|15fl . 
On the other hand, the non-inverse $ convention yields results for the activity 
coefficient given in Fig. l|16l) 

It is of interest to determine the activity coefficients by direct calculations be- 
cause the Kwong-Redlich and associated equations 42, n.29].|33] used to deter- 
mine activity coefficients for multicomponent mixtures are not easily extended 
to chemical reactions: part of the difficulty is that the component concentrations 



M ^2 



7~t-Ei,dC 



V + V(r 1 ,r 2 , ...Si....r M ). 
z — ' Zrrii 



(87) 
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(B) A, 



Figure 13: Single stage model with external potential interaction leading to 
chemical reaction, and which is similar to the elementary Bj0rn-Bjerrum theory 
of ionic reactions where (B) denotes the backward direction and (F) the forward 
direction of the Dimer reaction in eq.(0). 
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Figure 14: Double stage model with external potential interaction leading to 
chemical reaction, which goes beyond the single stage model and analogs where 
(B) denotes the backward direction and (F) the forward direction of the Dimer 
reaction in eq.JSJ). 



31 




Figure 15: Activity coefficients using the inverse $ convention due to boundary 
conditions. 




System concentration p 



Figure 16: $ convention determination of activity coefficients. 
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Figure 17: Availability coefficient ratios to illustrate the expected constancy of 
the availability coefficients of the excited species. 



in their thermodynamical mixtures are free to vary independently and may be 
fixed at any arbitrary value: this cannot obtain for reactive systems because 
of the relational dependence of the components via the equilibrium constant, 
and partial derivatives of thermodynamical quantities used to determine these 
coefficients demand free variation or fixing of concentration terms that is not 
realizable in a reactive system. Moreover there is the need to determine the 
variation of c(p) over the entire system density because of its theoretical sig- 
nificance linking $ and The next section provides a broad methodology, 
the exact methodology appears not to have been worked out but verification of 
the above would constitute another method of determining activity coefficients 
from rate measurements. 

3.6 Models based on simulation results 

The elementary rate constants have the form 03] 

fc 4 = A(T,r!)exp-^ (88) 

where E is termed the activation energy and Q are variables intrinsic to the 
reaction, such as impact and structure parameters. The object here is to relate 
the 7' reactivity coefficients to the 7 availability (activity) coefficient based on 
the first order single or two stage energy perturbation mechanism (Figs. $1'M 
Ir respectively. Refinements include perturbing the pathway itself, and intro- 
ducing a continuous potential field along the entire length of the reaction path- 
way. The object here is to present a broad framework where calculations to 



33 



1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





■ -o- y factor 


», 








\ 

\ 






\ 

Q. 

\ 

\ 

'•, 




\ 

\ 

X 




\ 




\ 




\ f 
\ 




o, 

\ 

\ 








o 



0.2 0.4 0.6 0.8 

P 



1.2 1.4 



Figure 18: Excited state activity for dimer before bonding for forward rate 
showing relatively lower change compared to the bulk activity coefficients (which 
would be the square of the values given in Fig. I|15|IV 
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any degree of accuracy might be attempted. For both these mechanisms, A\ 
and A* are the states of the reactants just prior to product formation. The 
first order perturbation lifts the degenaracy of the * levels relative to the vac- 
uum state; in the vacuum state, no singularities are observed in the potentials. 
Singularities due to the perturbation arises because the product and reactant 
states are distinct and distinguishable, and need not have the same activity in 
general; the form factor A(T, fi) is not altered to first order since they specify 
the type of reaction; only E is altered due to the external potential that modi- 
fies the ground state energies of the particles relative to the vacuum level. The 
availability coefficient is expressed as 

«ln7x (89) 

for any species X, and the excitation energy is written SB for species B.Foy 
what follows, e always refer to an energy term associated with the activity 
or availability coefficient. In the single stage model, the upper X* species is 
perturbed by an absolute energy amount given by 7x* in (|89(l : for the double 
stage model, the upper level is perturbed by the relative energy AE® FB ^ 2 , due 
to the singularity induced by the change of state from product to reactant or 
vice-versa. The lower state for both models are perturbed (lessened) by the 
same factor jx since this is the ground state availability . It will be shown 
that the 2-stage relative model is the more accurate and logical, based on a 
comparison with estimates of the activity coefficient for the monomer or atom 
A. The backward activation energy e b act for A2 — > A* 2 out — ► Product for either 
the 1 or 2 stage model can be written (superscript o refers to the vacuum state) 

-e'ct = -4 + *a 2 - 5A 2 , 0Ut (90) 

where we also write 5A2, ut — SAb.SAb. These terms refer to the upper state 
perturbation of the A^ out species. The forward model yields 

-<4t = -4 + £a + (-a- SA 2 , in (91) 

where 5A2_i n refers again to the upper energy state energy perturbations. Using 
to convert to availability, where SX = kT In ^ x , 190I91|I yield 

IaIa = • (92) 

lA 2 ,in 

On the other hand, the exact experimental condition in Observation^ can 
be written <f> = $' where the prime refers to the same expression in j' x where 

1 "fA 2 TA 2 (no\ 

1a 2 = = • (93) 

!A b lA b 7A 2 , ouf 



Hence, (|92I93|I leads to the exact result 

TA 2l< „ 7A 2tin 



7A 2 , out 7A b lA b 



(94) 



or SAb.SAb — bAi^in- In view of the fact that the product /reactant potential 
interface is completely different from the reactant/product potential interface 
for this hysteresis dimer, the result of l|94|) is indeed remarkable and may be 
stated as a kinetic principle: 
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Principle 1 The perturbed energy required to promote reactants to products 
at the reactant/ 'product potential interface is of the same magnitude as for the 
reverse transition of products to reactants at the product /reactant potential in- 
terface, even if these interfaces are discontinuous, i.e. are spatially and ener- 
getically distinct, for elementary reactions in equilibrium. 

Fig. (JJJ shows that for this study the interfaces are distinct;" 'reversible" ' path- 
ways have coincident interfaces. 

3.6.1 Two stage backward reaction model 
In view of (|94|l we can write the exact form 

i A = 05) 

a/7A 2 , out 

1M = (96) 

JA 2 , out 

If a complete characterization of the ja 2 species were known, then ja 2 out could 
be calculated from the theory presented in Sec. 1)3.3.1(1 especially expressions 
(I39|l . Some estimates or approximations may be made to determine the 7's. 
The reaction pathway for this 2 stage backward step model may be written 

Reactants{a) — > T.S.reactants(b) — > T.S. product s(c) (97) 

When considered in isolation, (normal vacuum state analysis) , there is continu- 
ity of the potentials between states (b) and (c) in H97|l . On the other hand, if the 
environment is considered where species are either reactants or products, singu- 
larities in the energetics would develop about the same arbitrary volume of the 
TS of the reaction coordinate; for instance at (b), two atoms in proximity A... A 
are separated by distance 5ta,A an d the mean activity of this state is strongly in- 
fluenced by the potential V(5ta,a) dependent on the 8ta,A distance, in addition 
to the surrounding non-participatory environmental molecules V (r others)' On 
the other hand, at the same vicinity (c) where the switch modifies the potentials, 
then the activity of the instantaneously formed molecule would be determined 
by V (r others), where V(5ta,a) now represents the potential of internal coordi- 
nates not connected to the activity. An apparent discontinuity arises according 
to this first order treatment, which has no apparent analog in the traditional or 
other QM theories computed at vacuum densities. For the two stage reaction, 
two types of approximations for activity coefficients may be made for the for- 
ward (F) 2A — ► A2 and backward (B) A2 — > 2A kinetic pathways and the results 
may be checked with the apparent activity coefficient values derived from the 
LJ single phase fluid. 

Case B: The first order perturbation energy E(^a 2 out ) is 

E(jA 2 ,out = 2e Ab - e* A2 (98) 

where Ab denotes the product atom just dissociated from the dimer at coordinate 
of species A\. The following approximation will be derived; 

SA 2t0Ut = ^^ (99) 

1A 2 



3G 



where jA b ~ JA- Clearly, SA2 tOU t serves as a '"retardant"' to the rate for posi- 
tive e's in the numerator of (|99l) . Fig. (|14|) can be correlated with (|97|l as follows: 
E B 1 is the activation energy to the level A* 2 , which consists of two degenerate 
states relative to the vacuum denoted TS(b) with molecule A 2 , and TS(c) which 
refers to the two atoms 2A after dissociation; in real time the molecule disin- 
tegrates along the trajectory (a)A 2 — > (b) — > (c) according to (JjJJJl. Unlike 
the atom, the dimer state may be characterized according to the inter particle 
A — A distance Ta-a, (as well as the magnitude of the external potential). The 
mean external forces acting on this dimer would lead to a potential which is 
to a first approximation relatively invariant, even if the internal potential and 
internuclear distance varies; thus we assign e(A2(a)) ~ e{A2{TS,b) . The po- 
tential energy is utilized to overcome the activation barrier; the boundary dC 
imposes a fixed density on the system which elevates the mean potential of the 
reactants relative to the vacuum ground state. During the transition to the 
state TS(b), external forces with a mean potential e(A2(TS) would operate on 
the dimer; r& = 1.20 for this reaction, where for a b.c.c. approximation for a LJ 
fluid at p — 0.70, the approximate nearest neighbor distance is 1.22, implying 
that the environment of A(TS) is very close to that of the bulk or generic A 
atom, i.e. a typical atom with activity coefficient 7^4, so that the transition 
state A 2 availability coefficient for the atom ja.ts can be equated with the 
bulk, ja.ts ~ 1a and the relative energy to be overcome from this TS state is 
AE° B 2 where AE° B 2 « (2e A {T.S.) - £a 2 (TS)) leading to 



exp 



A^, 2 
kT 



( J^J] . (100) 

\1A1AJ 



Because exp— AEb,i — exp~AE B 1 7a 2! we have the forms for the backward 
rate given as 



k B w k%j A 



1A 2 



lA,TSlA,TS 

_,/ 



k% lA ^ = k%^- y (101) 



But nnDHnUl implies 



' TA 2 _ 1 

7A, = — => lA 2 ,tn = jA 2 ,out ~ -T 

1A 2 , m $ 



leading to the useful approximation 

1A 2 « ^ (102) 

From JSEJ, ja 2 = lA 2 lA 2 ,out so that JHSl follows with j Ab ~ lA- Since 7A 2 ,in = 
lA 2 ,out, and i A i A = we get 

1A^% (103) 

Since the r.h.s. of H102I103() are available, we plot ja and -fA 2 in Fig. Q15[l. We 

expect 7^4 to be close or at least follow the trends of those derived from the Free 



37 



Energy estimates given in Figures (|10I12|> . especially Fig. (|1(J|> . The other two 
figures were plotted to show that the coefficient estimates are all greater than 
unity (> 1). As expected from the theory provided, the activity coefficients 
from the A res function show semi-quantitative agreement with the simulation 
results. The other A res functions which eliminate the dimer contribution shows 
markedly lower values; we would expect a slight lowering of value due to dimer 
interference, leading to the conclusion that simulation results with the above 
two models are plausible, and that A res is the more appropriate function to 
use for estimation. Since the adjacent atom in the dimer cannot contribute to 
the potential energy of the dimer contributing to its activity, one would expect 
its coefficient to be lower. The (B) reaction above has its counterpart in the 
forward (F) reaction, but because the transition state atomic activity coefficient 
is severely affected by its own force acting on the other target atom, it would 
not be reasonable to assume that the activity coefficient at that TS state due 
to the particle potentials can be equated with the bulk activity, especially when 
the internuclear distance at r/ = 0.85 corresponding to an energy of over 17.0 
LJ units! However it would be instructive to follow through the consequences 
of this approach so that comparisons with the activity coefficients derived from 
the literature may be made. 

F. Process Using the same arguments as for the B process, Fig. Ijl4|l (F process) 
gives the net elevation of the vacuum levels about AE"^ 2 such that 

6XP -— "felJ- (104) 

The ground state relative to the vacuum is elevated by e A per particle, so that 

E FA = E° FA - 2e A (105) 

where as usual the e's to the potential energy associated with the availability 
coefficients. The total activation energy AEpjot is then given by 

AE F ,tot = Ep \ + AE F2 (106) 

By the definition of the reactivity coefficients, 

~2e A -AE° F2 



IaI'a = CX P 



kT 



(107) 



and l(TU7|l leads to 



1 ' a7 ' a = Ja7a/^^\ (108) 

where the asterisked states are at the TS. Despite the arguments above, if we 
were to make the approximation 7^ w 7^4 (where in general , only a detailed 
perturbation theory would yield an appropriate value of 7^) then 7^ w j Az 
and 7* « j A yields 

« $ ( 109 ) 

IaTa 

Then ifTUSjl implies 

, , 1A1A , , 

lAlA = — (HO) 
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yielding 



1A « 7 ^V¥ (in) 

1A 2 « $7a 2 (H2) 

since here $ = $*. The above so-called $ convention plots for the activity 
coefficients is given in Fig. ()16|) where 7^, 7a 2 < 1, which contradicts the results 
of the free energy estimates given in Figs. (|1()I12(I . especially Fig. IjlOjl: this was 
to be expected from the assumptions made for state A*. Since &A% >OVj t can be 
determined from the more reliable (B) process, we might be able to derive an 
estimate for 7^7^ for the (A... A) TS state of two atoms about to dimerize where 
the (B) process assignment is used for the determination. We would expect 7^7^ 
to change relatively much more slowly compared to the ja and ja 2 variations 
over the system p if the activity coefficient is strongly influenced by the isolated 
large repulsive potential of the two atoms at the TS region. Evidence in this 
direction would serve as a prototype for describing ionic reactions along the 
classical ideas of Bronsted and Bjerrum. 



3.6.2 Single-stage considerations 

There is no reference state at the TS to compute energy differences. For in- 
stance, for the backward reaction 

1M = ^Ia * * (113) 

if 7^ w 7^4 as before for (B) (2 stage process). Then from $ = $' we get 7^1 
which is contrary to the simulation results. Similarly, for the (F) reaction, we 
have 7^ 7a = 7A , 7A and with the assignment 7^ sa ja 2 i we get i'a — ~~7g an d 

from l|96(l there results Y A2 w = 1 which is also not observed. Hence the 
absolute single-stage process are not expected although they seem to conform 
to the Bronsted-Bjerrum (B-B) form (where the z's are the associated charges) 

A z a + B z b _^ (A B ya+z<, products 

with k = k°A where A = 7 t j4Z " 7 - BZb with f referring to the upper level in 

the energy diagrams of Figs. (|13ll4|) or transition state of conventional theory. 
These observations suggest that the B-B form could well be described as an ap- 
proximation of another mechanism which does not have pre-equilibria transition 
states; an example being the two stage process below which subsumes the B-B 
equations. 

Ionic Reactions: The classical kinetic theory of salt effects applied to recent 
investigations by Sanchez et al. 03] gives the B-B equation for the rate 
constant k as 

k = k lAlB '" (114) 

7t 

where k$ is the reference rate constant for the solvent, 7^, 7s... are the reactant 
activity coefficients of reactants A,B,... and 71 is that of the transition state 
where if the charge of the reactants are z a ,Zb,... then the charge of the TS 
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denoted by f is z\ = z a + 2% + The activity coefficient 7,/(Debye-Huckcl 

limiting law approximation) for any species J conforms to 

log 7J = -Az 2 ! 1 ' 2 = -Qzj (115) 

where Q — AI 1 / 2 is a positive number, leading from (|1 14|) to the rate form 

log k 2 = log fcO + 2Az A z B I^ 2 (116) 

where a negative salt effect is expected for anion/cation bimolecular reactions, 
such as the recently studied reaction (with pz=pyrazine) 

k r 

[Fe(CN) 6 ] 3 ' + [R.u(NH 3 ) 5 pz] 2 + T f [Fe(CN) 6 ] 4 ~ + [Ru(NH 3 ) 5 pz] 3 + (117) 

A B CD 

where the forward reaction with rate constant kf shows a contradictory positive 
salt effect; other violations have been reported see ref. 4 of this citation]. 
The interpretation of these anomalies has been made using the theory of Marcus 
and Hush These workers introduced composite reactions, leading to a 

pseudo-elementary process such as 

fci 

A + B feTi PC 
PC ^ X^ sc 

SC -> C + D (118) 

for the overall reaction A+B — > C+D. Violations of the B-B formula l|114|) have 
been attributed to differences in activity coefficient properties of the precursor 
complex PC and activated state X'; SC is a postulated '"successor complex"', 
where the Markus-Hush ideas are also incorporated @B1 p. 15089]. It would be 
of interest therefore to frame theories and proposals for elementary reactions 
that might also explain some of the above results; this is attempted below. 



3.6.3 (i) Proposal of mechanism to explain positive deviation for 
forward reaction of ea. (|117|) 

Rewriting l|117(l with reactants C and D, products A and B with the charges, 
we have an elementary reaction 

C 1 - + D 3 + ^ A 3 - + B 2+ (119) 
The postulated reaction sequence route of an ideal model is 

C 4 ~...£ 3 + -> T.S 1 statc -> {^ 3 - + B 2 +} U State separated products P 

(120) 

In 1120(1 . C...D represents a pre- associated complex which retains the spectro- 
scopic details (e.g. in UV-IR range used for measuring concentrations) but 
which are close enough to ideally form an effective single charge of magnitude 
— 1. The T.S. is the lower (1) degenerate state relative to the vacuum, and P the 
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upper level u. In the presence of the solution dielectric, the first order pertur- 
bation according to the two stage model give a total pre-exponential total "/tot 
as 

1 ( 1A1B (u state) \ 
^ 7T .s.(l state) J 



Jtot = 7C...D 

for the overall rate k 2 = k\j t ot- Taking H115|l logarithms to base 10, results in 



log7tot = log7c...D + log7T.s. ~ log7A - log7s (122) 
= -1Q-1Q + 9Q + 4Q 

m +11Q (123) 

and so 

logfc 2 = logfc£ + llQ (124) 

which leads a positive deviation in contradiction to the standard B-B equation 
for cation/anion reactants. The above model had no intermediate forms; all 
charges were discrete etc.. AS a further refinement, one can write down inter- 
mediate forms of association with charge parameter —A and the I. state with 
charge — r leading to 

log k 2 = log k° 2 + Q(13 - A 2 - t 2 ) (125) 

It might be possible to write down intermediate forms of association along the 
lines of the above where the positive deviation is not a discrete multiple of Q 
above. 



3.6.4 (ii) Other mechanisms not in B-B form 

Large charged molecular fragments A Za and B Zb might not be large enough to 
be considered separate and discrete even at the TS region so that for the overall 
elementary reaction 

A z a + B z b ^ Products (126) 

with pathway A + B — > T.S. — > P, the two stage proposal yields the pre- 
exponential term j tot to hi where 

7m = iaib I (^r) • ( 12? ) 

If there is no further fragmentation, the charge on P would equal that of the 
T.S located at two centers, leading to jtot — JAIb with 

fc 2 = k° 2lAlB . (128) 

On the other hand, A and B if distinct would have the T.S. activity coefficient 
7t.s. ~ Ia-Ib and if P fragments, then the jtot multiplicative factor becomes 

Itot = {iaib) 2 Hp (129) 

where the jp is a product of the activity coefficients of the fragmented portions 
of the product states. 
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3.6.5 (iii) General harmonization of the B-B formulation with the 
present theory for free (non-associated) ions 

The ionic reaction with free ions 

#■ + B Zb -> Products (P) 
has the following reaction pathway for the two stage model 

A za + B z b ^ A Za ...B z " (T.S.(a)) AE=0 ^ acuum + B y-a+z b (T.S.{b)) -» P. 

. (130) 

AE act is defined to be the activation energy parameter from the free ionic state 
to the transition state a under vacuum conditions where there is bond formation 
and the smallest distance 5ta,b between two charge centers that are considered 
separate, and T.S. b is the smallest distance when they are not considered sep- 
arate and when a singularity of the external potential is applied as in the two 
stage model. Then the first order energy perturbation at the TS (between (a) 
and (b)) AE is 

AE = kTln( ^ A+B)Za+ '\ ) (131) 
and so the rate constant k% becomes outside of the vacuum state 

k 2 = fc 2 ° 1A1B .(7A.....7J'*) (132) 

where (7^4*0 ■■■jB z b ) is the TS species which has distinct charges relative to 
the dielectric medium and would appear as a product term of activity coefficients 
at the TS state, i.e. (7,4*0 ...7b*& ) = 7^4 ' -7 B ' • The energy of interaction 
of the A Za and B Zb would constitute a major portion of the activity coefficient 
of the pair, and would be expected to be relatively constant over a large range 
of system p compared to the bulk activity coefficients of the other species and 
reactants. Thus, under these conditions, we have 

fc 2 * fc°' 1A1B (133) 

!(A+By°-+ z b 

with fe;; = ^2 (jA*a ...jB z b ), which is the B-B equation for ionic reactions. 
Clearly, some demonstration of the relative constancy of ^A-a ■■■"fB z t> would be 
re-assuring. Obviously the potentials for the present dimer model and that of 
ionic reactions are different, but we can show a relative constancy for the analo- 
gous TS atomic pair, and so can expect a similar relative constancy to obtain for 
ionic reactions. We return to the previous system of ea. (|108|l where 7^ ■■■jB z b 
is equivalent to 7^7^ = {i*a---1a) f° r A = B. Then 

iWa = lAhl-il) (134) 

1 An 

A plot of this function is given in Fig. Ijl8|l . The analog of 7^ 9 is 7( J 4+s) s »+ z t • I n 

(|133|l . for A — B, jaJb = l\ would vary very dramatically, as is obvious from 

2 

Fig. H15(l . From p = to p = 1, varies by about 16, whereas (7x--7a) varies 
by 0.2. Thus, it is conceivable that the B-B expression may be derivable from a 
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strictly elementary reaction, based on the above estimates. Fig. (|17|l illustrates 
two plots, where the "Activity"' variable corresponds to using (|134f> to calculate 
7^4 with the assumption 7^ = ja 2 an d the '"Ratio"' variable is simply the 

reactivity and activity ratio using the inverse $ convention to determine 7,4- 

4 Conclusion 

The current form of the elementary reaction rate constant is still incomplete 
despite the depth and complexity of analysis valid for relatively low density sys- 
tems due to the utilization of traditional definitions and conventions. A more 
complete description would involve the '"reactivity"' coefficients and these coef- 
ficients are intimately related to the activity coefficients of the species involved 
in the reaction. Under stipulated conditions, such as obtain for the two stage 
model, explicit first order expressions may be derived connecting the reactivity 
and activity coefficients. A whole range of product and reactant states exists, 
and under an enumeration scheme incorporating the Gibbs equilibrium condi- 
tion, one particular state may be determined from another through a series of 
equations which also feature the bulk activity coefficients. The theory of ionic 
reactions as developed by Bronsted and Bjerrum may be subsumed by the sim- 
ple first order models provided here using a strictly elementary reactive process; 
in particular these models can also explain apparent violations of the standard 
Bronsted-Bjerrum theory involving composite reactions. The framework given 
here can be extended to other more complex particle reactions which are mod- 
erated by strong external fields, such as what might obtain in the plasma state. 
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